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Calculation of thermodynamic functions of the three-dimensional Ising ferro- 
magnet above and below critical temperature is performed in the approximation of 
sixfold basis distribution (p 6 model). Comparison with the results for the p 4 model 
indicates that dependence of the thermodynamic functions on the renormalization 
group parameter s becomes weaker. The optimal interval of the renormalization 
group parameter values is determined. 



Introduction 

Significant results in the description of the system thermodynamic properties 
in the vicinity of the transition point have been obtained by means of the col- 
lective variables (CV) approach. The method to deriving explicit expressions 
for the thermodynamic and correlation functions of the three-dimensional 
Ising model at temperatures both above and below critical temperature T c 
has been suggested within this approach. The calculations are performed 
with a non-Gaussian measure density. The measure is represented as an 
exponential function of the CV, the argument of which contains, along with 
the quadratic term, higher powers of the variable with the corresponding 
interaction constants. The simplest non-Gaussian measure density is the 
quartic one (p model) with the second and the fourth powers of the vari- 
able in the exponent. Then the sixfold measure goes containing the sixth 
power of the variable (p 6 model), and so forth. 

The results of the theory depend on the renormalization group (RG) 
parameter s due to an approximation of the Ising model partition function 
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calculation using the non-Gaussian measure densities. This dependence de- 
creases essentially if the non-Gaussian measure density becomes more com- 
plicated. Calculations of the correlation length critical exponent v within 
the p 2m models with m = 2, 3,4, 5 confirm this statement [1-3]. It has been 
established that the p 6 model provides an adequate description of the Ising 
model critical behaviour, in particular, the critical exponents, at the RG 
parameter values in the interval 2 < s < 4. 

Investigation of the p 6 model within the numerical realization of the 
CV method has been performed in [4]. Analytical derivation of the explicit 
expressions for the p 6 model thermodynamic functions is the subject of 
the present paper. The foundations for such kind of investigations have 
been developed in [5-9], where the quartic distribution was used as a basis 
measure. 



1 General relations 

The partition function of the three-dimensional Ising model within the six- 
fold measure density is given by 

\2n 



Z = 2 N [ exp[± ]T mk) P% P_ % + 2ni ]T u%9% + £ 3 (2 7 

J k<B k<B n=l 

x^- n 7ir E ■ ..^ k j, i+ ... + ,j{^f{dp)\ (i.i) 



where M 2 =l, M±=-2, M 6 =16, 4>(fe) = $(0)(1 - 2b 2 k 2 ), (3 = (kTy 1 is 
the inverse temperature, b is the effective interaction radius of the potential 
<3?(r) = Aexp(— r/b), $(0) = 8irA(b/c) 3 . Integrating (1.1) over p^ and 
ojj: with the indices B' <| k \< B (B = tt/c, c is the simple cubic lattice 
constant), we get an expression for the partition function of the p 6 model: 

Z = 2^e^'fexp[-^d'(%p_ r (1.2) 

J k<B' 

- e t^t^') 1 -' e - 2lP%i . . . P% j %i+ ... +% j{d P r . 

1=2 V '■ k u ...,k 2l <B> 

Here N' = Ns ~ 3 , s = B/B' = vrv^ft/c, 

d'(k) = a' 2 -(3$(k). (1.3) 
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Coefficients a 2 i depend on the ratio b/c and are given by the relations 

a' = lnQ(M), Q(M) = (l2sl) 1 /\~ 1 I (r ] '^'), 

a' 2 = {nsD^ftirf,?), (1.4) 

a 4 = 12s 3 C( V '^'), 

a 6 = (12slf/ 2 N( V ',a 

where the quantities rf = a/3sq^ 2 , £' = 8v 2 are the arguments, and special 

l5s Q 

functions C(?/,0 and 7V(r/,0 read 

C(r/,0 = -^4(^,0+3^(^,0, (1-5) 
^(r/,0 = ^e(»/,O-15^4(V,O^2(»/,O + 30^(»/,O- 

Here ^(V,0 = W.O = Jo°° ^e^~^dt. 

Using the method of layer-by-layer integration of the partition function 
in the phase space of CV, developed in [10], one can reduce (1.2) to the 
form: 

Z = 2 N 2^^ Z Q Z, . . . Z n (Q(P n )) N ^ J W f +l \p){dp) N ^ , (1.6) 

where N n = iV's" 3 ™, 

Z = [Q(M)Q(d)] N ' , Z 1 = [Q(P)Q(d 1 )] N \ ... , 

Z n = [Q(P n -i)Q(d n )] N », (1.7) 
i O) 

Q(Pn) = -(s 3 ^ — v) 1/4/ o(^,en), 

7T G{h n ,a n ) 
Q(d n ) = 2(24/4 n) ) 1 / 4 / (/ l „,a ri ). 

Hereafter, the arguments h n , a n are called basic: 

h n = d n (B n+1 ,B n )(6/a r i) 1/2 , a n = ^-a^ /(af ) 3 / 2 . (1.8) 
The effective measure density of the n-th phase layer Wq 1 ^ (p) has the form: 
wt\p) = exp[-- d n{k)pj:p^- (1.9) 



2 fe : l> 



1=2 ki,...,k 2 i<B n 



fci+-+fc 2 z J 
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Here B n = B's n . The intermediate variables r] n ,^ n are the functions of h n 
and a n : 

% = V6s 3 ^ 2 (h n ,a n )[C{h n ,a n )}-^ 2 , (1.10) 
in = ^s-V 2 N(h n ,a n )[C(h n ,a n )]-V 2 . 

The form of the special functions C(h n , a n ), N(h n , a n ) is being given by 
(1-5). 

Coefficients d n (B n+ i, B n ), are related to the coefficients of the 

n + 1-th layer by the recurrent relations (RR) [11-13]. The solutions of 
these relations [13] are used in the calculation of the system thermodynamic 
characteristics. 



2 Thermodynamic functions of the p 6 model in 
the regions of critical and limit Gaussian regi- 
mes (CR and LGR) above T c 

It is convenient to rewrite the model partition function as [14] 

Z = 2 N ZcrZlqr- (2.1) 
Let us consider Zqr given by 

Z C R = IIHTv ^) 1/4 /0(/in,«n)/0(r/n-l,en-l)] iV ". (2-2) 

n=0 17 U v7n-l)?n-lJ 

It should be mentioned that in (2.2) r/_i = 7/ , £_i = £' at n = 0. We repre- 
sent the right-hand side (RHS) of (2.2) in the form of an explicit dependence 
on the phase layer number n in order to calculate Zqr- 

In the CR region, the basic h n ,a n and intermediate r) n ,^ n arguments 
are close to their values at the fixed point. Therefore, functions of these 
arguments can be written as power series of deviations of basic arguments 
from their values at the fixed point (see [15,16]). Using the obtained repre- 
sentations for I (h n ,a n ), Io(f]n-i,£n-i), C^n^, £ n _i), we determine from 
(2.2) the partial free energy corresponding to the n-th phase layer: 

F n = -kTN n {f§l + MK-i ~ h^) + paK-i - « {0) ) + 
+^{K - fc(°>) + <^K " a {0) ) + <f'i(hn-i ~ h^) 2 + 
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+^K_! - aW) 2 + ^{h n - fc<°>) 2 + ^ - «(°)) 2 + 

V 5 (^-i - ^°))K_i - a(°>) + <Ps(K ~ h (0) )(« n - « {0) )}, 
/S = ln(^^)-ilnP 40 + ln/ *+ln/ **, (2.3) 
= 6 m + P 4m /4, m = l,2, 
^3 = -^2, = -^6, 

= C-^m-^4m/4 + ^L/8, <f! 3 = ft /2 - J? /2, 

^ = ^ 2 /2-^* 2 /2, ^ = 84-6x62 -P 43 /4 + V 41 V 42 /4, 

Expressions for the quantities occurring in y'i, <£>j are given in [15,16]. 

Hence, the partial free energy of the n-th phase layer F n is written as a 
power series of deviations of basic arguments from their fixed point values. 
The linear approximation for F n was used in [14]. In the present paper, 
as well as in the calculations within the p 4 model [5,6], the quadratics of 
the deviations are also taken into account. It allows one to compare the 
results of the calculations for p 4 and p 6 models. Let us note that quadratic 
terms of the RR do not contribute to the elements of the matrix of the RR 
linearization in the vicinity of the fixed point, and the eigenvalues E\ of this 
matrix and the critical exponent of the correlation length are the same as 
within the linear approximation for the RR. 

Let us find an explicit dependence of F n on the layer number n. Using 
the solutions of the RR, we get for h n and a n : 

h n = +c 1 H 1 ( U W)- 1 / 2 E? + c 2 H 2 { U W)- 1 E2 + 
+c 3j ff3(n (0) )- 3/2 ^3 n + c lC2 F4(u {0) )~ 3/2 £i n £ 2 n + 
+ ClC 2 tf 5 (n (0) )- 5/2 £ 1 n £ 2 2n + c 2 2 H 6 (u^)- 2 Ei n + 
+c\H 7 {u^)- l El n + clc 2 H 8 ( U W)- 2 El n E2 + 
+cl4H 9 (u^)- 3 Ef n Ei n , (2.4) 

a n = a ^ + c l L 1 (u^)- 1/2 El l + c 2 L 2 (u^)- 1 E^ + 

+c 3 L 3 (u^)- 3/2 E% + Cl c 2 L 4 (u (0) r 3/2 £M + 
+ Cl c 2 2 L 5 (u^r 5/2 E^E 2 2 n + c 2 L 6 (u^)- 2 E 2n + 
+c\L 7 {u^)- l E 2n + c 2 c 2 L 8 (uWy 2 E 2n E% + 
+c 2 c 2 L 9 (u(°Y 3 E 2n E 2 \ 
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where 

H3 = ^ _ , 

H4 = ^'»»W<?-f(l + wg'»<?), 

17._!^<» + ,«-^ fc m«g>), 

1 9 (0) o 

tf 6 = - ^g)), tf 7 = ^(!*<«M? - ^6), 

_ V6 (o) SgW^ ^6 (o) 3a<°> 

Ll = "is ^ 31 — ' L2 = T^ 32 " — ' 

Ls " Ts — ' (2 - 5) 

1/15 (o) \/6 (o)x 

wfiV 5 (o) (o) (o)x 
L 7 = y(j« tu, «2i " — '), 

Considering (2.4), we rewrite the partial energy of the n-th phase layer as 

Fn = -^ 8 ^[/S+/S(^v i ^+/S('' (0) r 1 ^+ 
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+/S(« (0) )- 3/2 C3^ + f^ R (u^r 3/2 c lC2 E^E^ + 

+/S (^r^dc^r + /S(u(°))- 2 ci^ + (2.6) 

+/S(^°))- 1 c^i 2n + /S(u(°))- 2 c 2 c 2 ^^ + 
+/S(n (0) )- 3 c 2 c^ 1 2 ^ 2 2 ™]. 

Here, 

fcR = ^(^3 + <Pi/Em) +L m ((p4 + (p 2 /E m ),m = 1,2,3, 
/g> = H 4 (<p 3 + ip 1 /(E 1 E 2 ))+L 4 (<p 4 + <p 2 /(E 1 E 2 )) + 

+2H 1 H 2 (<p 3 + tp'jiE^)) + 2LiL 2 (v? 4 + ip' 2 / (EiE 2 )) + 

+(#iL 2 + ff 2 £i)(P6 + (p's/iEiEh)), 
/g> = H 5 (ip 3 + Vl /(E 1 E 2 ))+L 5 ^ 4 + ^ 2 /(E 1 E 2 )) + 

+2(H 1 H 6 + H 2 H 4 )(ip' 3 + ^/(^i^ 2 )) + 2(LiL 6 + L 2 L 4 ) x 

x (994 + v 2 l(E 4 E 2 2 )) + (fTiL 6 + fl" 6 Li + # 2 L 4 + #4^2) x 

x^ + ^/^E 2 )), 
f§ R = H 6 (<p 3 + cp 1 /E 2 ) + L 6 (<p 4 + <p 2 /E 2 ) + H 2 (<p' 3 + ip[ /El) + 

+L 2 {^ 4 + if' 2 /E 2 ) + H 2 L 2 (<p' 6 + (p' 5 /E 2 ), (2.7) 
/g> = H 7 (ip 3 + i Pl /E 2 ) + L 7 (ip 4 + ip 2 /E 2 ) + H 2 (ip' 3 + ip' 1 /E 2 ) + 

+L 2 (<p' 4 + <p' 2 /E 2 ) + HiLi(ip 6 + if' 5 /E 2 ), 
f§l = H 8 (cp 3 + ip 1 /(E 2 E 2 )) + L 8 (ip 4 + i P2 /(E 2 E 2 )) + 2(H 1 H 4 + 

+H 2 H 7 ){ip' 3 + tp'jiEfE^) + 2(LiL 4 + L 2 L 7 )(ip' 4 + 

+(p' 2 /(E 2 E 2 )) + (F1L4 + F 4 Li + #2^7 + #7^2) x 

x(^ + ^/(E 2 E 2 )), 
fg ] R = H 9 (cp 3 + ip 1 /(E 1 E 2 ) 2 ) + L 9 (i P4 + cp 2 /(E 1 E 2 ) 2 ) + (2H 1 H 5 + 

+2H 2 H 8 + Hi + 2H 6 H 7 )(<p' 3 + ( p' 1 /(E 1 E 2 ) 2 ) + (2L 4 L 5 + 

+2L 2 L 8 + L 4 + 2L 6 L 7 )&' 4 + V 2 l{E 4 E 2 f) + (fTiLg + 

+H5L1 + H 2 L 8 + H 8 L 2 + H 4 L 4 + HqL 7 + H 7 Lq) x 

x(^; +¥>;/(£! £ 2 ) 2 ). 

The quantities n^,it;i ^ were determined in [13]. Let us note that in the 
expressions for /i n and a n we can neglect a qualitatively new term propor- 
tional to E 3 which arises within the p 6 model considered (£3 is not essential 
as compared to E\ or E 2 ). 
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Summing up expressions for F n (2.6) over layers of the phase space from 



n = to n = m T , we obtain 



F CR - F + F C R, 
Fq = -kTN[\nQ(M)+\nQ(d)], (2.8) 

t CR - CR ' CR + ^CR T ' 



where 

,(o) ,(1) -1/2- i-, 

F<°> = -fcTiV^- 3 [-^_ + Wo cir^i 
CR L l-s~ 3 

, fcR^o 1 C2oE 2 f§ ] R <pQ 3/2 c 30 E 3 f^ ) R <Po 3/2 c 1 Tc 2 oE 1 E 2 

1-E 2 s- 3 l-E 3 s~ 3 l-E^s- 3 

, fSW'^ircl^El fS^o 2 4 El 
1-E^s- 3 1-E^s- 3 

fSW3f3, , &o 2 cjT 2 c 20 E!E 2 
+ l-^s- 3 + 1 - E%E2S- 3 

f^cir^EfEl 

,(0) ,(1) -1/2 , 

F« = kTN's- 3 "* [-^ + /cf?(/?0 1° + 

CR L l-S- 3 1-^iS- 3 

I fcFt^O 1 /o 1 

,(2) 4. /C 3 ^0 3/2C 30^- +1 , 

^ - -kiNs [ - i _ E ^_ 3 + + 

f (4) -3/2 f p m r +l ,(5) -5/2 , o p 2(m T +l) 
I -/CR^O jQC20-fc 2 -/aR^O J0 C 20^2 , 

, /c 6 RV 9 Q 2c 20- g 2 (mT+1) , fcR^o' ' fo c 2$E™ T+1 

I -E^s- 3 l-E 2 1 E 2 s~ 3 

f (9) -3 f 2_2 F 2(m T +l) 

I JCR^O J0 c 2Q rj 2 i 

+ 1 - £ 2 £ 2 2 s -3 J - 

To calculate these expressions, we have used the formulas: 

s -3(m T +l) = T 3^-3m 0) ^ jq) 

d = rci/9$(0), c 2 = C20 (/?$(0)) 2 , c 3 = C30 (/?$(0)) 3 , 
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where mo, ci, c 2 o, C30 and /o, <po were defined in [14]. The dependence of too 
on s for the p 6 model under consideration is plotted in figure 1 (solid curve) . 
Here the dashed curve corresponds to the p 4 model. It is easy to notice that, 
when s = s* (s* is the value at which h n turns to zero at the fixed point; 
s* =3.5862 for the p 4 model and s* =2.7349 for the p 6 model), the values of 
toq for p 4 and p 6 models coincide. 




Figure 1: Dependence of the quantity mo on the RG parameter s for the p 6 
model (solid curve) and the p 4 model (dashed curve). 

(2) 

Further, we put Fq R = Oat r < 1, since E 2 < 1, -E3 < 1, and m T is large. 
Taking E2 , E3 into account gives rise to terms characterising corrections to 
scaling. As a result, the free energy of the CR region takes the form: 

F CR = -A;TiV'[ 7o + 5o-7f R)+ T 3 1, 

,-3r fcl 1 fSW'hrE, ff R ^ 2 c 20 E 2 

l-E 3 s~ 3 l-E 1 E 2 s~ ;i 

A 5 ).- 5 / 2 - „2 _ P P 2 f(6) -2 2 p2 /(7) -1-2 2p2 

1 /cr^o cic 20 tEiE 2 JcrPq c 20 -fc 2 , Jcr^o c l r ^1 , 

f(8) -2~2„ ^2^2^ j?(9) -3~2_2 „2p2p2 
, JCR^O C 1 C 2QT £g^2 Jcflyp CiC 2 qT E 1 E 2 l 

+ l-EfE 2S -z + 1-^5-3 J ' 1 j 

-(0) 

<f = lnQ(M) + lnQ(d), 7 f R)+ = = (^-f, 



/o 



AO) ,(1) -1/2, ,(7) 1,2 

+ _ JCR 1 JCR^O JO , JCR^O JO 

1 1-a- 3 l-^iS-3 ^l-^s-3- 
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Note that 7o,<5o are the functions of temperature, since they are ex- 
pressed in terms of c\,C2o,csq and Q(d). Let us extract the temperature 
dependence in these quantities. 

Near T c we have for c\ 

~(o) , ~(i) 
c\ = c\' + c\ 't, 

i (0) -1/2 

~(o) rr n ^ , (o) 1/2 . (o) , _i . a 4 ^i2 

q = vi[i-/o + «l2Vo + 4jYwo + ^^0^2 + 

^JV] (212) 

' (0) -1/2 „ / (0) -1 
~(1) y r Q 4^1 2 Vo , 3a 6^1 3 Vo 1 

1 1L (PM0)) 2 (PMoyf J ' 

for C20 and C30 we get, respectively, 

(0) , (1) , (2) 2 

C20 = c 2( / +c k 20 't + c^ 't , 

' (0) 1/2 

(0) r (0)/-, , s 1/2 (0), -1/2 «2 V 21 ^0 

c 20 = ^2[-W)-»2l(l-/0W -V23^0V0 + ^ ^ + 

' ' (0) -1/2 

+ " 4 +Wgo j (213) 
(/3 C $(0))2 (/? c $(0))3 

' (0) 1/2 ' o ' (0) -1/2 

(1) y r Q 2^21 ^0 , 2 «4 , g^ggVp 1 

20 /3 C *(0) (/3 C ^(0)) 2 (&*(0))3 J ' 

' o ' (0) -1/2 

J2) _ T/[ Q 4 , 3q 6^23^0 i. 

c 20 — K 2l , I ~ „ T 



(&$(0)) 2 (/3 C *(0))3 



(0) , (1) , (2) 2 

C30 = c 30 ; + c 3C /r + c 30 ; r , 



,(o), 



(0) T/ r , 3/2 (0) n (0) 2 W 31 (/3 . 

3 30 = ^[-^ " V?0 U 32 -W - /o)«31 + ^ |^ + 

' (0) 1/2 ' 
+ M^L_ + ^^ ]; (2 . 14) 



-30 



' (0) ' (0) 1/3 9 ' 

.(1) = y^ Q 2^3lVo | 2q 4^32V ] 3q 6 



/3 C $(0) (/? c $(0))2 (/3 C $(0))£ 

,(0), 1/2 



.(2) T/ r Q 4 ^32 , 3a 6 

'30 — v 3[ , „ ~ „ -r - 



(/3 C $(0)) 2 (&<K0)) S 
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The values of /3 C <3?(0), the correlation length critical exponent v = lns/lnEi, 
exponents of the corrections to scaling Ai = —In E% /In Ei, A2 = 
— In E3 / In Ei , and coefficients of the expressions for c~i , C20 , C30 are given in 
tables 1,2. In the present paper, the numerical calculations are performed 
at b/c = 1 and arithmetically averaged Fourier transform of the potential. 

Table 1: Values of /3 C 4>(0), u, A 1: A 2 , , c ( p for different s. 



s 


&*(0) 




Ai 


A 2 


g(0) 


~(i) 


2 


1.1204 


0.619 


0.653 


5.061 


0.7602 


0.0109 


2.5 


1.1405 


0.634 


0.552 


3.963 


0.7637 


0.0103 


2.7349 


1.1506 


0.637 


0.525 


3.647 


0.7641 


0.0100 


3 


1.1628 


0.640 


0.503 


3.379 


0.7630 


0.0096 


3.5 


1.1882 


0.645 


0.476 


3.038 


0.7570 


0.0089 


3.5862 


1.1929 


0.645 


0.473 


2.994 


0.7555 


0.0088 


4 


1.2165 


0.648 


0.460 


2.821 


0.7471 


0.0082 


Table 2: Coefficients in 


equations for C20 (2.13) and C30 (2.14). 


s 


c (0) 

L 20 


c (1) 


c {2) 

c 20 


c (0) 
c 30 


c {1) 

c 30 


c (2) 

c 30 


2 


-0.2302 


-0.0949 


0.0075 


0.2401 


0.0762 


-0.0161 


2.5 


-0.3297 


-0.0893 


0.0080 


0.3266 


0.0698 


-0.0147 


2.7349 


-0.3820 


-0.0872 


0.0083 


0.3642 


0.0654 


-0.0140 


3 


-0.4449 


-0.0851 


0.0085 


0.4056 


0.0603 


-0.0133 


3.5 


-0.5737 


-0.0813 


0.0089 


0.4821 


0.0516 


-0.0120 


3.5862 


-0.5972 


-0.0807 


0.0089 


0.4952 


0.0503 


-0.0118 


4 


-0.7142 


-0.0776 


0.0090 


0.5581 


0.0443 


-0.0108 
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Expressions for the quantities Vi, V2, V3, ipo, vff occurring in (2.12) 
(2.14) are presented in [13,14]. The coefficient 70 can be written as 

(0) , (1) , (2) 2 

7o = To +To t + To T > 

AO) f (2) -1 (0) p f (3) -3/2 (0)p 

_ „-3r jCR , /Cfi^O c 20^2 Jcr^O C 3Q £3 
TO " S l?3^3 + l_^ 2S -3 + l_£ 3S -3 + 



+ 



1 - Els" 3 J 



,(1) -l/2-(0) p ,(2) 1 (1)„ ,(3) -3/2 (1)„ 

..(1) „-3r /gflgg C l £;1 1 JCR^O c 20^2 JcR^O c 30 ^3 



^fcR^O 2c 20^ c 20^2 1 



+ ^ T _^ % (2.15) 

(2) -3, fe^l , fSW^l 1 f { cW 2 ^ 

,(4) -3/2, -(0) (1) -(1) (0) w p 
, /cfl^O ( c l c 20 + c l c 20 
1-^^25-3 

f (5) -5/2, (0) (0) (1) (0)x 2 ^p j?2 

, /Cfl^O ( 2c l c 20 c 20 + c l ( c 20 j )hlE 2 

+ 1 - EiS 2 2 S -3 + 

+ /g^o 2 ((4 1 o ) )^+24°o ) 4 2 o ) )^ , f&Vo\f?l% , 



1 - £|s-3 1 - E\ 



. /g^0 2 (ci 0) ) 2 4^^2 , /g^0 3 (4° ) ) 2 (^ ) ) 2 ^ 1 
1-^25-3 ^ \-E\E\s-* J ' 

For 5o we obtain 

4 0) = In Q(M) + In Q(d,T c ), 

6 (i) = _^ ={1 _ q) ^ {0):F2{h ^ a): (2 . 16) 



4 2) = -^(l-qf((3 c ^)) 2 [^(h Cl a)-T A {h c ,a)] + 
a 4 
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/fi 



a/6 r i n ~ , . , _.. a/6 a' 6 _ 1 + s~ 



2 



,|a 2 -A*(0)(l -,-)], «-^^, • 

Coefficients of the expressions for 70 (2.15), 5q (2.16) at different values of 
the RG parameter s are given in table 3. 

Table 3: Values of coefficients in equations for 70 (2.15) and 5q (2.16). 



~M ~( 2 ) a( 2 ) 

7o 7o 7o °0 °0 °0 



2 


0.1205 


-0.3160 


-2.5291 


0.2711 


-0.3574 


0.4810 


2.5 


0.0757 


-0.2343 


-2.9320 


0.3237 


-0.4510 


0.6464 


2.7349 


0.0624 


-0.2062 


-2.9801 


0.3423 


-0.4862 


0.7128 


3 


0.0510 


-0.1798 


-2.9798 


0.3606 


-0.5222 


0.7827 


3.5 


0.0362 


-0.1415 


-2.9314 


0.3905 


-0.5836 


0.9075 


3.5862 


0.0343 


-0.1360 


-2.9236 


0.3953 


-0.5938 


0.9287 


4 


0.0267 


-0.1135 


-2.9039 


0.4175 


-0.6420 


1.0319 



Hence, the free energy of the CR region reads 

F CR = -kTN'[^ CR) +7 ir + 7 2T 2 + 7 f R )+ r 3H, 
lf R) = 7f (2-17) 
7! = 

7 2 = 7f+e 

The numerical values of coefficients 7q C '" R ' ) , 71 , 72 an d quantities 7 + , c v = 
(c^V/o)^ occurring in are given in table 4. 

Knowledge of Fcr allows one to calculate other thermodynamic func- 
tions of the system in the CR region at T > T c . For the entropy Scr, 
internal energy Ucr and specific heat Cqr we get 

s CR = fcAry°)(^) + CoT + 4^)V-«], 
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Table 4: Coefficients ^q CR \ 71, 72 and quantities <7,,7 + ,7 contained in 
^(CR)± different va i ues f the parameter s. 



s 


(CR) 

To 


7i 


72 




7 + 


7 


2 


0.3915 


-0.6734 


-2.0482 


1.4412 


-0.3170 


0.7382 


2.5 


0.3994 


-0.6852 


-2.2856 


1.2722 


-0.7757 


0.5244 


2.7349 


0.4047 


-0.6924 


-2.2672 


1.2097 


-0.9831 


0.4188 


3 


0.4116 


-0.7020 


-2.1971 


1.1462 


-1.2229 


0.2899 


3.5 


0.4267 


-0.7251 


-2.0239 


1.0414 


-1.7450 


-0.0277 


3.5862 


0.4295 


-0.7298 


-1.9949 


1.0250 


-1.8496 


-0.0973 


4 


0.4442 


-0.7555 


-1.8720 


0.9517 


-2.4345 


-0.5148 



U C r = fcT^[7i+«ir + 4 Cfl) V- a ], 

C C r = kN'[cv + cf R)+ T- a ], (2.18) 

smcR) = 7 (^) +7l , C0 = 2( 7l + 7 2 ), 4 CS)+ = -*nlF*> + , 

Ul = 2 72 + 7l , cf R)+ = -3i/(3y - lhf R)+ , a = 2- Sis. 

In the region of LGR, the expression for the part Z^gr of the partition 
function (2.1) reads 

Zlgr = j "exp{-i [dm T ( k ) ~ dm T ( B m T+1 ,B mT )]pj:p_j:- 

k<B mT+ i 

- 2™ £ «sps-^ E rlT^F-i- (2-i9) 

fc<B mT + i k<B mr+ i 



«=2 k U ...,k 2 l<B mT +l 

h 1+ -.. + kJ(dp) NmT+1 (^) NmT+1 - 



To calculate Zlgr it is convenient to select two regions of the wave vector 
values [14]. The first, transition region, corresponds to the values of k close 
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to B mT ; the second, Gaussian region, corresponds to small values of the 
wave vector (k — > 0). Hence, we have 



ZlGR - Z LGR Z LGR- 



(2.20) 



The contribution to free energy from the phase space layers following 
the point of exit from the CR region is 



7 (i) 
LGR 



-kTN f T RT 3u , 



-3m 



fLGRi(m), 



(2.21) 



Itr = cIJtr, Jtr = ^2 s 

2 1 1 

fLGRi(m) = In - + - In 24 - - lnC(r] mT+m , 6n T +m) + 

7T 4 4 

+ lnJ (/i mr+m+ i,Q; mT+m+ i) +hi/ (??m T 



where rao is the nearest integer to rh . The quantities rh , rj mT+m , £ mT + m , 
hm T +m+i, a m T +m+i were determined in [14]. The plots of rh (s) and anal- 
ogous dependence mo {s) (p 4 model) are represented in figure 2. 



IB 







\;. m o 











4 5 

S 



Figure 2: Behaviour of rh (p 6 model) and m (p 4 model) with the increase 
of the parameter s. 

Introducing an infinitely small external magnetic field h = hb'H (l^B is 
the Bohr magneton), we can write the part of the free energy corresponding 

to Z \jQfl ^S 

F% R = -kTN'f'r^ - (3Nlth 2 T- 2v , f' = 4f\ 
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fLGR 2 = -\\n2A + \ + \\nu mTl _ l - l -HG+^)- (2.22) 
- ^^(V-i'^'-i) " Gs 2 + (G S 2 ) 3 / 2 arctan[(G S 2 )- 1 /2] ) 
7 4 + = c; 2 7 4 + /(/3$(0)), 7 4 + = ^°/(2G). 
Here, 

m' T = m T + rh + 2, u m ^_ 1 = ^^(/^(O)) -2 , 
G = (u^m^iMh^^,^^)]- 1 ~ q. 

A general expression describing the contribution of long-wave fluctua- 
tions to the free energy (LGR region) reads 

Flgr = -kTN'f LG RT 3v -(3Nl4 + h 2 T- 2v , (2.23) 
fLGR = clfiGR, J LGR = fTR + I ■ 

The values of Jtr, /', 74" are given in table 5. 

Table 5: Values of Jtr, /', 74" ■ 



s 


fTR 


/' x 10 5 


it 


2 


0.6529 


0.4749 


2.7055 


2.5 


0.8142 


0.0102 


3.0870 


2.7349 


0.8824 


0.2155 


2.1737 


3 


0.9541 


0.0659 


2.1841 


3.5 


1.0756 


0.0091 


2.2058 


3.5862 


1.0950 


0.0067 


2.2098 


4 


1.1822 


0.0016 


2.2307 



The entropy, internal energy, specific heat corresponding to the LGR 
region are defined by the relations 
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Slgr = kN'u^r 1 -*, 

U LGR = kTN'4 LGR) T 1 - a , (2.24) 
Clgr = kN'4 GR r- a , 

(LGR) n 

H = 3vf LGR , 
c (lgr) = 3 ^^ 3z/ _ ^ f LGR 



Contributions to the thermodynamic functi- 
ons of the model from the critical and inverse 
Gaussian regime (CR and IGR) regions below 



C 



The free energy at T < T c can be written as [6,10] 

F = F + F CR + F IGR , (3.1) 

where Fq = —kTNln2 is the free energy of the system of N non-interacting 
spins, Fqr is the contribution to the free energy from the short-wave fluc- 
tuation phases of the spin moment density (CR region), and Fjqr is the 
contribution from the long- wave phases of the fluctuations (IGR region). 

The number /x T of the CV phase space layer, separating the short-wave 
and long-wave phases of fluctuations, is an important characteristic of the 
system. It is determined from the equation 



r (o) 



-(0) 



5. (3.2) 



Here corresponds to the fixed point of the RR, r^+i is determined from 
the solutions of the RR equations (see, for example, [13,14]), 5 is a constant 
(S < 1). In the present paper we put 5 = 1 (see [8]). Let us write an 
equation for fj, T 

I T I c 1 E^ +1 = f , (3.3) 

the solution of which is 

In I r I ln/o - lnc^ 0) 

^ = - W + W-l, Mo= lnEi ■ (3.4) 
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We need to sum the partial free energies over the layers of the CV phase 
space to calculate Fcr. Extracting an explicit dependence on the layer 
number, using relations (3.3) and 

s -3(A*r+l) = | T |3^ S -3 M0) S -3 M = c 3 ; (3.5) 

we obtain 

F CR = -kTN'[^ CR) - 71 | r | +72I r | 2 - 7 f r I 3 "]. (3.6) 
Coefficients 7q ,71,72 are determined in (2.17), 

tF^ = cfr - , (3.7) 

f (0) f (l) -1/2- ,(7) l, 2 

7 I-S- 3 l-^S- 3 

The value of 7~ is given in table 4. 

The entropy, internal energy and specific heat of the system correspon- 
ding to the CR region read 

S CR = kN'[sWV - co | T | +uf R) -\ r M, 

C/ Ci? = fcTAr'[ 7l - «i | r | +uf R) ~\ r | 1_a ], (3.8) 

C C i ? = fciV'[ C o-cf R) -|rr a ], 

^3 =3^73 

4 Cfl) ~ = 3i/(3i/ - 1)7^" ■ 

Let us calculate now the contribution to the free energy from the IGR 
region 

F IGR = -kTN's'^+V in[v^Q(p^)] _ fcTlnZ^+i, (3.9) 

where 

4. s 3„(Mr) 

^ Q(P " J = ( ^ra» 1/4/ »("-^»- 

^ r +i=/exp{-- ^ ^+i(^)PfcP_fe- 
fc<B MT+ i 

-E^<ti E pfe -p^ + ... + fe,](dp) JV - +1 - 

«=2 k 1 ,...,k 2l <B„ T+1 



18 



Consider the first term on the RHS of (3.9). Making use of the relations 
r> r = -7> T /3$(0), f^ T = Ul + E^), 
u^ T = u^(/?<i>(0)) 2 , =^o-rf4°i ) ^r 1 1 (3-10) 

w = ^(^(o)) 3 , = - /o^i^r 1 , 



w 



we find 



>/6 



7 /Try — ' /ir V" «V 



U = ^|^ 3/2 iV(V,a MT )[C7(V,^)]" 3/2 - 
The first term on the RHS of (3.9) is equal to 

s-^ +1 hn[V2Q(P, T )] = lg \ r l^ + ^+Din^^, (3 . 12) 

7 9 = 4%, % = H( , r t s h UflT n J V %(?K^J]- 

To find the second term on the RHS of (3.9), we need to calculate Z^ T+ ±. 
The coefficients occurring in it equal 

<W+i(*0 = 7> T+ iS- 2 ^ +1 ) + # 2 , g = 2/?<l(0)6 2 , 

V+i = -^+1/3^(0), f^+i = 2/ , (3.13) 
u Mr+ i = u Mt+ i(/3$(0)) 2 , u^ t+1 = ip - foip 1/2 w ( 2i, 

W/Xr + 1 = ^ T + l(/?^(0)) 3 , W^ + l = V0 - /O^OW^. 

Let us perform the change of variables 

p s = p' n + VN<*>6 i (3.14) 

in the expression for Z^ T+ \ in order to extract the free energy related to the 
ordering that has appeared in the system. Here < a > is being determined 
from the extremum condition [17,10] 

<"> 2 = io "niV '- i+ ^- < 3 - i5) 

a 6 
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b, = 



i 



1 + 



6a^ +1) | d^ +1 (0) 
5(a^ +1) ) 2 



Simultaneously, we include in the treatment a constant external field h = 
Hb7~L, which sustains the separated average moment, and separate from the 
sums over k the terms with k = 0. We obtain 



Z^ T+1 = exp(-/3F (T + (3F h ) J dp exp{(3VNhp - ^ T+ i(0)p; 



l(Mt+1) 



■>r+l) 



!>T+1) 



a, 



(i"r+l) 



fc<B, 



Mt+1 

AT,,. 



+P3/£»5+P4Po}(dp) V 



+1- 



Pol* 



(3.16) 



The prime on the sum over A; means that k 7^ 0, 



= H I ^ T+l(0 ) I ^l—l^+xC-l + 62) - 



(Mr+l) 



25 (a 



(Mt+1)\3 



18 ( a (^+ 1 ))2 



^ T +l(-l + ^) 2 



(3.17) 



PF h = 



^^F^yi ( _ 1+fe2)1/2- 



We have for the integrand coefficients of (3.16) 



<W+i(*0 = 4 I <V+i(0) 



10 {a { t +l) ) 2 

3 „(Mr+l) 



(-1 + 62) + qk 2 



= cl I r | 2 ' /?$(0)[4r V+ i - H!k±l(_i + 6 2 )] + # 2 



3 w 



(3.18) 



■>t+1) 



10 (Mr + 1) / (Mr+1) , 

1 (-1 + 6 2 )^[l + + 6 2) ], 



&^ +1 >=a? T+1) [l + 5(-l + & 2 )], 
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The quantities pi are being determined by the equations 

3 X 4 

6 (^r+l) 

x E ^fci ' ' ' / ? fcA 1 +...+fc 4 ~ 5!jv L^/tvttT x 

a { ^ +1) 

x E ^ ■ ■ ■ ^ 5 ^feiH hfes fi!/V 2 

- E Pfc! ■ ■ ■ Pk 6 feiH h*^' 

fei<S M r+l ^ T+1 ^<-B MT+ l 

Pl ~~2j^~; ^ p k p ~k~ 3 \ N ~ 2^ Pk 1 Pk 2 p k 3 x 

~^nv2 E ^fei •••/ 5 fc 5 % 1+ ...+fc 5 ' ( 3 - 19 ) 

V+i fel <i? MT+1 

X E PkiPk 2 Pkz^ki+k 2 +k3 ~ ARM 2 ^ ^1 " P kJ*ki+-+kV 

E ^Vfes 



X 

ki<By, T+1 

>r+i) 



P4 = ~^T7i k £ nP -"' 

Mt + i k<B^ T+1 



u 2 

The dependence of the quantity Bq = 4f Mr+ i — ^ ,/ T+ * (~1 + ^2) on the 
parameter s is plotted in figure 3. The dashed line shows the dependence of 
the analogous quantity 4/o on s for the p 4 model. One can see from figure 
3 the quantity Bq ocurring in d^ T+ i(k) (3.18) (therefore, and d^ T+ \{k)) is 
positive for all s. 
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1 2 3 4 5 6 7 8 

S 



Figure 3: Dependence of the quantities Bq = 4f^ T+ i — X w^+i + ^ 2 ) 
model) and 4/o (p 4 model) on s. 



Expanding exp{po +P1P0 + P2P0 +^3/°o +P4Po) m ser i es an d restricting 
ourselves to the terms of second order, we integrate (3.16) over p^ with 
k 7^ using the Gaussian basis distribution. Gathering up the series over 
the averages with respect to the Gaussian distribution in the exponential, 
we obtain [15] 

Z^ T+1 = exp(-PF a + PF h -PF m ) []' ( t \,J 1/2 x 

h<B^ +1 °W+iW 

x J exp(A lPo + i 2 Po + A 3 pl + Mpl + A h p\ + igpg + 

+A 7 pl + igp^dpo, (3.20) 



where 



1 1 / „(Mt+1) h (p T +i) 

-PF m = ^ T+1 [4(e +1) + L Il)Il2 + i_ l3X 

X (^ + 2?1 3) + (a ^ +1))2 (^ + ^(^ + ^)) + 
v 5 1 ; 8 V 45 2 V 4 3 

I.(^T+l)„(/ir+l) T -7-2-7- 

-^1(^7 + — ^— ;J) 



4 v 3 2 

ix = /3v^V/ i + i< /2 +1 [-J 1 (^ +1) + ^f +1) ^i) + 
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+ 2 ( ¥ + ~ )} + 4 ( 15 + 



;; 1 (i 1 4 



"2 | <Wi(0) | +j|^±__L(_i + b 2 ) + ^[-J!(6i^ +1) + 
a 6 

+^^x 1 ) + + (^ +1) ) 2 (§ + ^) + 

+ 2 ( ¥ + Jl( Y + — )] + 4 ( 15 + 

+ 1 + ^))+4^ +1) 4^ +1) (f + 

+6 (, T+ i) a (, T+ i) ( | + Ji( | + (321) 



(X T + 1 



+ I 1( f + ^I 2 )) + ^ — x 2 ], 

_± r 1 /7>r + l) , ^6 T \ , H T , 

+ + + X 2 + 

a (^+i) 6 (^+D (> T +i) )2 j 4 r 3 7^^^ 
+ 9 J3 + L V^ ( f + Jl( f + 16 JlJ2))] ' 

fJ, r -\-l 

b (v-r+i) ( Mt+ i) (Mr+ i) X 3 XiX 2 

H 2 — 2l J +ft 5 a 6 ly + — JJ> 
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~ i A!" +1> , i 4"* +1 ' (>>'i" +v ' f- 

„(Mr+l) „(^r+l) 
+ ^X 2 (^ +1) + ^X 1 )], 

I>t + 1),>t + 1) 
7 _ °5 Q 6 T 

288< /2 +1 

8 23047V3 t+1 2 - 

Expressions for X m occurring in (3.21) can be found via the transition 
to the spherical Brillouin zone and integration over k G (0, -B^+i] [17,18]. 
We have 



J i = w— E' i4tm=^I^ + i(o)' 10 



fc <t; r+1 T 3 

(-„(m t +1)n 2 



X 



or 



s 2(//r+l) 

/3$(0) 



a x . (3.22) 



Here, 



£ = 3 Xr ~ ar f anXr , x r = , \ (3.23) 



5 V+ 1 4r W -f|^(-l + 6 2 ) 



The quantity Z 2 can be represented in the form: 

J 2 = V # 2 (r) = a 2 , (3.24) 

V (/3^(0)) 2 

where 

1 ^/ 6s 2 ^ T+1 ^ 
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2 

1 (-1 + b 2 ) + l}- 1 [sin(B ftr+ ir) - B^ T+1 rcos(B^ +1 r)}, 



3 V+i 

a 2 = a? + 6e?(l + ei), (3.25) 

6 

e i = ^ . 

20 



e2 



— [sin(7r\/2) — 7r\/2 cos(7rV2)] ~ 0.034861. 
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The other 2j(Z = 3, 4, 5, 6) are determined by analogous relations 

s 2Z(m t +1 

WW) 



Ti = Y.9\r)= (3-26) 



with ^ 

^ = «i + 6ei(l + ^n). (3.27) 
The values of a m {m = 1, 2, 3, 4, 5, 6) are given in table 6. 

Table 6: Numerical values of a m . 



s 




«2 


"3 


«4 


a 5 


Q 6 


2 


0.3456 


0.1893 


0.0488 


0.0151 


0.0050 


0.0017 


2.5 


0.2981 


0.1406 


0.0313 


0.0083 


0.0024 


0.0007 


2.7349 


0.2821 


0.1258 


0.0265 


0.0067 


0.0018 


0.0005 


3 


0.2664 


0.1122 


0.0223 


0.0053 


0.0014 


0.0004 


3.5 


0.2412 


0.0918 


0.0165 


0.0036 


0.0008 


0.0002 


3.5862 


0.2373 


0.0888 


0.0158 


0.0033 


0.0008 


0.0002 


4 


0.2201 


0.0764 


0.0126 


0.0025 


0.0005 


0.0001 



Let us perform the change of the variable po m (3.20) 

Po = Po ~ < a >, (3.28) 

which cancels terms proportional to odd powers of po in the exponent of the 
integrand. We obtain 

= exp(-/3F; T+1 ) J exp[(3VN Po h + Bp\ - *Lp* - 
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■§2Po]d PO , (3.29) 



5 5 

a lf- , \ , ^/l+l r «4 . 2 

+ a V + [y + «i«2(- + 

5 , ... 5 _ ,_ w u +i . 

+-(-1 + b 2 ))\ - -r llT+1 a 1 a 2 [u ltT+1 + «i] + 

^r+l r«6 of 04 o|«2 

32 L 45 2 1 3 4 )l 



16 3 2 



-o^ [I" 1 + 6 2)(«i - n^ T+1 a 2 (- )) - 

3«V+i 2 3 w^ T+ i 

-(7 - 56 2 )r V+ i« 2 ] - \ M^8y) " 5 MU + 

+4rV + i-^^(-l + ^))vr- 1 ] + i-^}, 
3 to^+i 3 3 

^ | r \ 2v [3<l>(0)f^ +1 B, (3.30) 

| r |" (/34»(0)) 2 S [jG, 

1 - 2^(«^ + i + — «0 + ^ T^" 1 + ^ >< 
. V+i > , 1 /«i«2 . a 3 5 

X( T + ^ + 2 { — + T )] " 2 a2( ^ +1 + 

+ ^^ } + 8^ [ 15 + ai( ^~ + 2 ( T + ^~ ))] + 



H — - h ai , 

2f MT+1 L 12 u 8 3 ,b 



+ + 5f MT+1 « 2) ] - + §(-1 + 
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+ai{—a 1 a 2 + y)J, 

- + i 1 U^ + l wl T + 1 a 3 OLXOL2. 

The quantities B, G, D are given in table 7. 

Table 7: Values of B, G,D,<a >(°),r+,r- for different s. 



s 


5 


G 


D 


< a >(°) 


r+ 


r- 


2 


0.8296 


0.0153 


0.0002 


3.0942 


2.6052 


0.2970 


2.5 


0.8289 


0.0206 


0.0003 


3.0031 


3.8147 


0.3121 


2.7349 


0.8277 


0.0231 


0.0003 


2.9669 


2.9709 


0.3211 


3 


0.8267 


0.0260 


0.0003 


2.9309 


3.3248 


0.3318 


3.5 


0.8266 


0.0316 


0.0003 


2.8730 


4.0679 


0.3528 


3.5862 


0.8268 


0.0326 


0.0003 


2.8640 


4.2071 


0.3565 


4 


0.8282 


0.0376 


0.0004 


2.8238 


4.9254 


0.3748 



Using the method of the steepest descent for Z^ T+ i (3.29), we find 



Z^r+i = y^ry exp[-/?F; r+1 - NE (p)}. (3.31) 

Here p is the extremum point of the expression 

Eo(p) = Dp 6 + Gp 4 - Bp 2 - (3hp (3.32) 
arising in the exponential of the integrand of (3.29) at the substitution 

p = s/Np. (3.33) 

For Eq(p) at h = we obtain 

E (p) = -s~ 3 ^) s ~3 Eo , 



F 2& ( 1 l i/l l 3 f ^ BD ) i f »r + iBG 



2 L 3f^+iBD, 



^- 1+ 3V 1+ 2^g^)- 
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Having (3.12) and (3.31), we can calculate the contribution to the system 
free energy at T < T c from the long-wave phases of the spin moment density 
fluctuations (3.9): 

F IGR = -kTN\ 1 { t ) + 1 r > )\rf\ 

l ( t ] = 7 ff + 7p = 473 ( ^ ) , i^=% + %, 

s- - 5_ . 5_ n"? 
1p = c„7 P , lp = 2 r Mr+H«i + -r^+ia 2 ) - — x 

x(^ + i + + + a\a 2 {\ + + 6 2 ))] - 

-\f, T+1 a ia2 [u, T+l + ^±l ai ] + + | x (3.35) 

x( T + ~ jJ + 16 [ T + ~ J 3^ x 



x [(-1 + b 2 )( ai - V+i«2(^ - ^^^)) - (7 - 56 2 ) x 

2 6 w^ T+ i 



! 1, r 1 + 4 ^-ffe(- 1 + ^) 1 , 1 £ 
xrv +1 a 2 ] - - ln[ ] + - - -, 

-,<cr> „Z- <cr> - <<T> 771 

73 - c */73 , 73 - 

The quantity 73 1 " determines the free energy after the exit from the CR 
region, and "y^ a> determines the free energy of the ordering. The values of 

lg-> Ipi 7s* T \ lz a> are given in table 8. 

The entropy, internal energy and specific heat of the system correspon- 
ding to the IGR region read 

SlGR = Sfj, T + S <a> , UiGR = Ufj, T +U <a> , Cigr = Cfj, T + C <a> , 
S v = -kN' \T\ l - a v%\ U n = -kTN' \r\ 1 ~ a u { ^\ (3.36) 
C v = kN'c^ ] I r \~ a , u^ = 3u^\ 4")= 3i/(3i/ -1)7^. 

The index 57 can take two values: fi T < a >. 

4 The order parameter of the system 

The mean spin moment is the order parameter of the model investigated. 
It is related to the presence of the nonzero value of the CV po at which the 
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Table 8: Values of 7 9 , 7p, 73 , 73 



s 


% 


% 


1^ 


7 3 <CT> 


2 


-0.3024 


1.0386 


0.7362 


1.7618 


2.5 


-0.0869 


1.0269 


0.9399 


2.0456 


2.7349 


0.0039 


1.0227 


1.0265 


2.1572 


3 


0.0986 


1.0179 


1.1164 


2.2808 


3.5 


0.2579 


1.0076 


1.2655 


2.5154 


3.5862 


0.2832 


1.0056 


1.2888 


2.5563 


4 


0.3967 


0.9955 


1.3922 


2.7538 



integrand of (3.29) has an extremum. Performing the substitution (3.33) in 
that integrand, we obtain 

ZfiT+l = e-PKr+iy/N J e- NE ^dp, (4.1) 

where Eq(p) is given in (3.32). Owing to the factor TV in the exponent in 
(4.1), the integrand has a sharp maximum at the point p corresponding to 
the equilibrium value of the order parameter. The value of p can be found 
from the extremum condition 9E Qp P ^ = or 

6Dp 5 + 4Gp 3 - 2Bp - j3h = 0. (4.2) 

In the case h = we obtain a biquadratic equation, which is reduced by 
means of substitution 

p 2 = y (4.3) 

into the equation 

6Dy 2 + 4Gy - 2B = 0. (4.4) 

Extracting the temperature dependence, we obtain the equation for the 
mean spin moment < a >= p = ^fy: 

<a > = \rf<a >(°), /? = i//2, 
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< a > 



(0) 



c y 2 (/?<i(o)rvv /2 <^> (0) , (4.5) 



<^> (0) = [^(-i + f + p^)] 172 - 

The value of < a >(°) is given in table 7. 

The susceptibility per particle \ can be found from equation (4.2) by 
differentiating it with respect to H and using the relation x = P-B 9 g^ > '■ 

= _ (4g) 

30Dp 4 + YlGp 1 - IB 

Separating the temperature dependence in the coefficients D, G, B (see 
(3.30)), we obtain a final expression for the susceptibility. 

5 Thermodynamics of the system in the vicinity 
of the phase transition point 

Having calculated the contributions to the system free energy from the short- 
wave and long-wave modes of the spin moment density oscillations both 
above and below T c , we can compute the total free energy 



F = 



F + F CR + F LGR , T > T c , 
F + F CR + F IGR , T < T c , 



the entropy, internal energy and specific heat. 

We obtain the total free energy of the system at h = taking (2.17), 
(2.23) and (3.6), (3.35) into account: 



-£:TiV'[7o + 7iT + 72T 2 + 73 + T 3! 1, T > T c , 
-kTN'[ 10 + 7l | r | +72 | t | 2 +73" | r \ 3 % T < T c , 



where (see table 9) 

7o = 7f R) + ^ln2, 

7 3 + = -7i° R)+ + hGR, (5.2) 

(C'R)- . (ht) , <<r> 

73 = -73 + 7/Gfl, HGR = 73 +73 ■ 

The coefficients 7^ can be written as a product of the universal part 7^ (ta- 
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Table 9: Coefficients 70, 7 3 and 73 . 



7o 7 3 7 3 7 3 7 3 



2 61.1798 
2.5 61.1878 
2.7349 61.1930 

3 61.1999 
3.5 61.2150 
3.5862 61.2179 

4 61.2325 



0.9699 2.9033 

1.5898 3.2734 

1.8654 3.3020 

2.1770 3.2783 

2.8206 3.1856 

2.9445 3.1704 

3.6167 3.1179 



1.7599 5.2680 

2.4612 5.0675 

2.7650 4.8944 

3.1073 4.6793 

3.8086 4.3013 

3.9423 4.2448 

4.6608 4.0180 



ble 9) and non-universal factor (? v depending on the microscopic parameters 
of the Hamiltonian (the lattice constant c, the effective interaction radius b 
and the value $(0) of the interaction potential Fourier transform at k = 0): 



7 3 ± = 



7 3 + = "7 + + Ilgr = -7+ + /+ / , (5-3) 



:,<a> 



73 = -7 + Ig + lp + 73 ' 
We have for the entropy, internal energy and specific heat of the system 

f kN' + c r + r 1 ""] , T > T c , 

\ kN [ S W -co \ t\-Us \ t I 1 ""], T<T C , 

{ kTN'l^+mT + u+T 1 -*], T>T C , 

\ kTN'[ ll -u 1 \r\-u^ \r\ 1 ~% T<T c , { ' } 



C = 



kN'[co + 4T- a ], T>T C , 
kN'[c + Cz I r \-% T<T C 
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with the coefficients given by the relations 



u 



,(0) 

± _ 



7o + 7i, 



3„-;± 



c„u 



uf = 3vyf, 



± 



u u 3 > 



ui = 272 + 71, 



(5.5) 



+ 



c 3 c^ 

C,.Co 



-3 - ^3 > ^3 - 3l/ ( 3 ^ - 1)7^- 

The formula for the specific heat can also be rewritten as 



C/kN 



a 



A ± = c 3 acf, 



r \- a +B ± , 
B ± = c . 



(5.6) 



The plus and minus signs correspond to T > T c and T < T c , respectively. 

The system susceptibility per particle at infinitely small values of the 
external field H at T > T c ( X = - ^$$1* , see (2.23)) and T < T c (see 
(4.6)) is given by 



X 



1 T »(o)' J > Jc ' 



-7 V^B 



$(0)' 



T<T C . 



(5.7) 



Here (see table 7), 



v 74 



; 2 { 



10 G 2 



3 D V 2 5 V 2 CP 



7 



= 2v. 



-1 



(5.8) 



Plots of the temperature dependence of the system free energy F/N, 
entropy S/kN, specific heat C/kN, mean spin moment < a > (4.5), sus- 
ceptibility x (5-7) at s = 2,2.5,3 are shown in figures 4-8. 

The corresponding plots of the thermodynamic characteristics calculated 
within the p A model, with the confluent corrections [9] being taken into ac- 
count, are also given here (dashed curves). Comparison of these plots shows 
that the dependence of the thermodynamic functions on the parameter s for 
the p 6 model is weaker than for the p 4 model. The dependence of F/N on 
s for these two models is represented in figure 9. 
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Figure 4: The temperature dependence of the free energy of the system for 
different values of the RG parameter s within the frames of the p 6 model 
(solid lines). For the comparison we show the free energy of the system in 
the quartic basis distribution approximation with allowance for confluent 
corrections [9] (dashed lines). (1) s=2, (2) s=2.5, (3) s=3. 
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Figure 5: Dependence of the system entropy on r (Notations are the same 
as in figure 4). 
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Figure 6: Specific heat of the system (Notations are the same as in figure 
4). 
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Figure 7: Mean spin moment < a > at T < T c (Notations are the same as 
in figure 4). 
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Figure 8: The temperature dependence of the susceptibility x (Notations 
are the same as in figure 4). 
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Figure 9: Behaviour of F/N as a function of s for p 4 and p 6 models. 



35 



Let us note that the calculations performed are best suited for the in- 
termediate values of s, close to the quantity s*, at which h n turns to zero 
at the fixed point (s*=2.7349 for the p 6 model and s*=3.5862 for the p 4 
model). The use of the difference form of the RR based on a non-Gaussian 
measure density works especially well for this region of s. At small values 
of s (s — > 1), some complications arise when the unit element is extracted. 
In this limit, the RR should be represented as the perturbation series with 
respect to the Gaussian distribution (h n is large at s — > 1, and expansions in 
h~ 2 , a n can be used [19,20]). There also exists an upper limit for s. At large 
s, one must take into account the correction due to the potential averaging 
[10, 21], which increases with s. 

The point s ~ s* corresponds to the beginning of the u(s) curve stabili- 
sation [3]. Table 10 contains the values of the critical exponents u, a, (3, 7, the 
exponent of the scaling correction Ai = — In E2/ In E\, and the ratios of the 
critical amplitudes A + / 1 A~ ,T + /T~ and their combinations P = ^[1 — ^], 
R+ = A + T+/[sl(< a >(°)) 2 ] at s = s* for p 4 and p 6 models. These values 
are in agreement with the data obtained within the field theoryapproach 
(FTA) [22-27] and high-temperature series (HTS) [28-32]. 

Conclusions 

The method for the calculation of the three-dimensional Ising model ther- 
modynamics is developed in the sixfold distribution approximation. Both 
temperature regions above and below the critical value of T c are considered. 
The main distinguishing feature of the approach is the separate allowance 
for the contributions from the short- and long-wave fluctuation phases of 
the spin moment density to the free energy of the system near T c . Within 
the framework of the p 6 model, we obtained the explicit expressions for the 
critical amplitudes of the thermodynamic functions of the three-dimensional 
Ising ferromagnet and calculated the coefficients of the free energy, the uni- 
versal characteristics (the critical exponents, the ratios of the critical am- 
plitudes). Calculation of the free energy, entropy, specific heat, mean spin 
moment, susceptibility was performed for different values of s. Comparison 
of the results obtained within p 4 and models indicates that the depen- 
dence of thermodynamic functions on the RG parameter s is weaker for the 
p % model. 

The values of s close to s* are optimal for the method presented. Obtain- 
ing analytical expressions for critical amplitudes and system thermodynamic 
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Table 10: Values of the critical exponents, ratios of the critical amplitudes 
and their combinations for s = s* (s*=2.7349 for the p 6 model and s*=3.5862 
for the p 4 model) obtained by means of the CV method. Data calculated 
within the field theory approach (FTA) [22-27] and high-temperature series 
(HTS) [28-32]. 



Quantity 


Model 


Ref. 


data 


P 4 


P b 


FTA 


HTS 


V 


0.605 


0.637 


0.630 




0.638 


a 


0.185 


0.088 


0.110 




0.125 


P 


0.303 


0.319 


0.325 




0.312 


7 


1.210 


1.275 


1.241 




1.250 


Ai 


0.463 


0.525 


0.498 




0.50 


A+/A- 


0.435 


0.675 


0.54, 0.48 




0.51 


r+/r- 


6.967 


9.253 


4.77, 5.12 




5.07 


p 


3.054 


3.711 


3.90, 4.03, 4.2, 


4.72 






0.098 


0.162 


0.059, 0.052 




0.059 



characteristics as functions of the Hamiltonian microscopic parameters is the 
advantage of the method developed. The leading critical amplitudes for the 
specific heat and other thermodynamic characteristics are represented as a 
product of the universal part, independent of microscopic parameters, and 
the non-universal factor, which depends on these parameters. 
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